A METHOD AND APPARATUS FOR DETECTING AND LOCATING NOISE 
SOURCES NOT CORRELATED 

FIELD OF THE INVENTION 
The present invention relates to detecting and 
5 locating sources of noise in the general sense, using 

sensors that are appropriate for the nature of the noise 
source . 

The invention relates to a method of detecting and 
locating noise sources disposed in a space of one, two, 
10 or three dimensions and optionally correlated with one 
another, and presenting emission spectra of narrow or 
broad band. 

The invention finds particularly advantageous 
applications in the field of locating sources of noise 

15 optionally accompanied by echo and coming, for example, 
from vehicles, ships, aircraft, or firearms. 

BACKGROUND OF THE INVENTION 
In numerous applications, a need arises to be able 
to locate in relatively accurate manner a source of noise 

20 in order to take measures to neutralize it. Numerous 
solutions are known in the prior art for acoustically 
locating noise sources. The main known solutions make 
use of techniques for correlating signals delivered by 
detection sensors . 

25 Those techniques present the drawback of being 

particularly sensitive to interfering noise occurring in 
the environment of the measurement sensors. Furthermore, 
it must be considered that those techniques constitute 
specific methods that are adapted to each application 

3 0 under consideration. 

The technique in most widespread use involves 
antennas having a large number of sensors (several 
hundred) and a large computer system implementing beam 
forming so as to aim in a given direction in order to 

35 increase the signal-to-noise ratio. That method does not 
make any a priori assumption concerning the number of 
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sources and any possible correlation between them, which 
leads to a loss of resolution. 

OBJECTS AND SUMMARY OF THE INVENTION 
There therefore exists a need to have a general 
5 method of detecting and locating noise sources in space, 
when the number of noise sources is small and is known or 
overestimated . 

The invention seeks to satisfy this need by 
proposing a method of detecting and locating noise 
10 sources by means of sensors adapted to the nature of the 
noise source, the method presenting low implementation 
costs . 

To achieve this object, the method of the invention 
consists : 

15 • in taking the time-varying electrical signals 

delivered by the sensors each signal 3^(t) delivered 

by a sensor being the sum of the signals Sj emitted by the 
noise sources (X^) ; 

• in amplifying and filtering the time-varying 
20 electrical signals as taken; 

• in digitizing the electrical signals; 

• in calculating the functional 

fin-^, Tij, n^) = X-^/ci 

with the coefficients -R^^.^ being a function of the vectors 
25 rij. giving the directions of the noise sources; and 

- in minimizing the functional in such a manner as 
to determine the directions of the noise sources. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Various other characteristics appear from the 
3 0 description given below with reference to the 

accompanying drawing which shows embodiments and 
implementations of the invention as non- limiting 
examples . 

Figure 1 is a diagram showing the principle of the 
35 detection method of the invention. 



Figure 2 is a diagram showing a detail 
characteristic to the method of the invention. 

Figure 3 is a diagram showing the method of locating 
two noise sources using two sensors. 
5 MORE DETAILED DESCRIPTION 

As can be seen in Figure 1, the method of the 
invention consists in locating noise sources X^, X2, .../ 
Xj, . . . , X^ where j varies over the range 1 to M, the 
sources being distributed in space and each emitting a 
10 respective signal Sj with j varying in the range 1 to M. 
The method of the invention consists in locating the 
noise sources Xj using sound wave or vibration sensors Y^^, 
¥3, . . . / Yi, . . . , Yn where i varies over the range 1 to N, 
each delivering a respective time-varying electrical 
15 signal s^, S2, - • - / s^, . . . , s^* 

The method consists in taking the time -varying 
electrical signals s^(t) delivered by each of the sensors 
and representative of the sums of the signals emitted 
by the noise sources X^ . The signals 3^(t) received on 
2 0 the N sensors on the basis of the sum of the 

contributions of the various sources is written as 
follows : 

where i = 1 to N, r^^ is the distance between the noise 

25 source X^ and the sensor Y^, and c is the speed of sound 

in the ambient medium. 

The term A^^ represents the attenuation due to 

propagation together with the sensitivity factor of the 

sensors and is expressed as follows: 
30 A,^ = B,C(r,^) 

where i = 1 to N and j = 1 to M, where is the 
sensitivity coefficient of sensor and where C(r^^) is 
the attenuation coefficient due to propagation over a 
distance . 
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The sensors are associated with respective 
electronic units (not shown) for amplifying and lowpass 
filtering the signals they pick up. The sensors are 
preferably matched in modulus and phase so that their 
5 sensitivities are identical. Thus, = G for i = 1 to N. 

Advantageously, in order to facilitate implementing 

the antenna of sensors as defined above, the sensors 

are placed relatively close to one another. 

Consequently, for remote sources, the distance r^j is of 

10 the order of the distance , i.e. the distance between 

the center of gravity of the sensors and the source X^. . 

Thus, attenuation becomes a function of the distance r^. 

only with C(r^^) = C(r^) , with i = 1 to i\r and j = 1 to M. 

It can be deduced therefrom that: 
15 A^j. = G.C(r^) = a(rj) 

where i = 1 to N and j = 1 to M and: 



where i = 1 to N. 

Since the amplitudes of the sources X^- are unknown, 
20 the following equation can be written as follows, 
integrating the term a.(r^) in S^i 



where i = 1 to N. 

Using Fourier transforms, the expression for the 
25 signals s^(t) becomes: 

(1) Si (CD) = 2:Sj(a>) .e ^ 

whe re i = 1 t o N 

where s and S are the Fourier transforms of s and S 
respectively and where O) is angular frequency. 
30 This first equation (1) relates the received signals 

to the distance ^±j, i.e. to the positions of the sources 



As can be seen in Figure 2, other relationships can 

be expressed, associated with geometrical considerations 

enabling the distances r^j to be related to the unit 

vector n^, which determines the direction defined by the 

center of gravity of the sensors and the source 

generating the signal S^. The position of the sensors is 

defined by the vector constructed from the positions of 

the sensors and the position of their center of 

gravity. A development restricted to the first order of 

r^j then provides : 

(2) r-. ^ - <n^., > 

where i = 1 to 2^7" and j = 1 to M, and where . > is the 

scalar product . 

Thus, by replacing r^^ by the approximate expression 
given in (2) and integrating the phase term: 

e 

which depends only on the source X^. in the magnitude 
S ^{(O) , equation (1) can be written: 

(3) Si (CO) = j^s^m 

where i = 1 to N. 

This relationship can also be expressed in matrix 
and vector form: 

(4) Si(Q)) = i;s^((o) .rj.(co) 

with, for ith coordinate of the vector T^i 

(TjOi = e <^ 

where i = 1 to 
Or indeed: 

(5) s(CO) = T.sm 
where T = matrix having the general teirm: 
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When the sources are not correlated, the signals 
Sj can be determined from the signals of the vectors 
zij. Cross-correlation functions between and for i ^ 
j are then minimized. 
5 Once the minimization operation has been performed, 

after determining the directions n^, it is also possible 
to discover the magnitudes S^. 

If N = M, i.e. if there are as many sensors as 
sources, then the system (5) can in general be inverted, 
10 If N > M, the problem can be reduced to a square 

system by premultiplying by: 

i.e. by the conjugate transposed matrix of T. System (5) 
then becomes: 
15 ^T*.8{(0) = ^T* .T.SicO) 

I.e. 

(6) Sico) = (""T* .D'^.^T* .8{Q)) 

With the signals S expressed formally in this way, 
the correlation coefficients R^j between the sources i and 
20 J are calculated formally by: 

where can also be calculated formally from frequency 
magnitudes , giving : 

jZ\Sii(0)\\\s,{(x)fd(0 



(8) J?,,. 



j!:|s,(co)fdco.j::|s,((o)|'do) 



25 The function for minimizing is then: 

(9) firij^, n^., n^) = 

where the coefficients Rj^^ are functions of the vectors 



When the signals S^, Sj are received with comparable 
30 amplitudes, the denominators of R^^j are of the same order 
of magnitude and can then be replaced by 1 without 
spoiling the positions of the minimas. Calculating the 



r^j. can then advantageously be performed in the time 
domain, when the range of variation in possible delays is 
small . 

Once the directions defined by the vectors have 
5 been determined, it is also possible to find the 

magnitudes from equation (6) . Such a technique thus 
makes it possible to determine the natures of the sources 

The description below with reference to Figure 3 
10 gives an example of detecting and locating two noise 

sources X-^, X2 that are not correlated (M = 2) , using two 
sensors Y^, (N = 2) . 

In the frequency domain, the electrical signals s^^, 
S2 delivered respectively by the sensors Y^ and Y2 and 
15 representative of the sum of the signals S-^, emitted by 
the noise sources X^^, are expressed as follows: 



where 



T.. = ^ 



2 0 is the propagation delay of the signal emitted by source 
i prior to reaching sensor j . 

Inverting this system leads to: 
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The cross -correlation function between the source 
2 5 signals and S2 is written: 

^12^ = Jt: 5^(co) .S^{co)e''^ dco 

for the delay T. 

Replacing 5i((0) and ^2(05), it becomes: 
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whence 
N(Q)) = 



ri2(t) = J 



\D«i)f 



dm 
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and 



|D(a))|^ = 4sin^— (x^j^ + 
2 
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A sample (but sub-optimal) solution in this case 
consists in optimizing the numerator N only. 

The cross-correlation then be approximated 

by: 

1^12 (X) = j!» iV(CD)dco 

Replacing N((D) by its value an expression is 
obtained which is a function only of the y^j corresponding 
to the autocorrelations and cross-correlations between 
the measured signals and s^i 

ri2(T) ^ -Yll('T + ^12-T^22)-Y22(^+'«^ll -'^2l) + Yl2('«^ + T^ll-'^22) + Y2l('^ + '^12-'^2l) 

It is recalled that the distance r^^ can be 
approximated by: 



r • . = r • 



Thus, replacing jt^^. in this approximate expression 
and integrating the phase term 
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in SjiO)) finally leads to an expression of the estimator 



of which is as follows: 



r,2(T) = - Y 
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where zij = the unit vector of (OX^) with i = 1, 2. 
However : 

<nT,c-,> = - — cos 01 
2 

D 

< 223^,02 > = — cos 9^ 
D 

D 

<n2,C2> = cos 02 

where the distance between sensors is written is D. 
Then: 

ri2('0 - -Yiil^-t + ^ cos 02 j 

-y22{j + ^ cos 01 j 

+ Yi2|t^ + 



+ (cos 01 + cos 02 )j 

+ Y21^'^ + (cos 01 + cos 02)^ 



The functional to be minimized relative to (0i, 02) 
is thus: 

10 - ri\(x)dT 



sign ambiguity between 0^ and 02 is removed by- 
analyzing the half -plane containing the sources and 
assumed to be known a priori . 
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The invention is not limited to the examples 
described and shown, since various modifications can be 
made thereto without going beyond this ambit. 



Without ftirther elaboration, it is believed that one skilled in the art can, using the 
preceding description, utilize the present invention to its fullest extent. The preceding 
preferred specific embodiments are, therefore, to be construed as merely illustrative, and not 
limitative of the remainder of the disclosure in any way whatsoever. Also, any preceding 
examples can be repeated with similar success by substituting the generically or specifically 
described reactants and/or operating conditions of this invention for those used in such 
examples. 

Throughout the specification and claims, all temperatures are set forth uncorrected in 
degrees Celsius and, all parts and percentages are by weight, unless otherwise indicated. 

The entire disclosure of all applications, patents and publications, cited herein are 
incorporated by reference herein. 

From the foregoing description, one skilled in the art can easUy ascertain the essential 
characteristics of this invention and, without departing from the spirit and scope thereof, can 
make various changes and modifications of the invention to adapt it to various usages and 
conditions. 



